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We develop a comprehensive theory of phase for finite-dimensional quantum systems. The only 
physical requirement we impose is that phase is complementary to amplitude. To implement this 
complementarity we use the notion of mutually unbiased bases, which exist for dimensions that 
are powers of a prime. For a d-dimensional system (qudit) we explicitly construct d + 1 classes of 
maximally commuting operators, each one consisting of d~ 1 operators. One of this class consists of 
diagonal operators that represent amplitudes (or inversions). By the finite Fourier transform, it is 
mapped onto ladder operators that can be appropriately interpreted as phase variables. We discuss 
the examples of qubits and qutrits, and show how these results generalize previous approaches. 
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I. INTRODUCTION 

The standard formalism of quantum optics is usually 
presented in the context of the harmonic oscillator, where 
both position and momentum are represented by un- 
bounded operators with eigenvalues in the real numbers. 
Systems living in a finite-dimensional Hilbert space were 
studied originally by Weyl and also by Schwinger 0, 
but except for some relevant exceptions (for a complete 
review see 0), they have received the attention they 
rightly deserve only after becoming one essential ingredi- 
ent in the development of the emerging field of quantum 
information [j, |5| • Indeed, the promise of futuristic tech- 
nologies like safe cryptography and new "super comput- 
ers" , capable of handling otherwise untractable problems, 
relies on the ability to control the quantum states of a 
small number of qubits 0, Q . 

In the modern parlance of quantum information the 
concept of phase for a d-dimensional system (or qudit) 
is ubiquitous. However, in spite of being a primitive of 
the theory, this notion is rather imprecise and, roughly 
speaking, three quite distinct conceptions can be dis- 
cerned. 

In the first, phase is considered as a parameter and 
the problem is reduced to the optimal estimation of the 
value of the phase shift undergone by the qudit under 
certain operations |8( . Although very operational in style, 
it accommodates perfectly the practical requirements of 
typical applications handled in quantum information. 

In the second, a semiclassical approach is adopted, and 
the phase is assumed to be linked to the geometry of the 
state space: for example, for a qubit this space is the 
well-known Poincare sphere and the phase is identified 
with the angle that a state representative makes with 
the Z axis [9j. This pictorial understanding of phase as 



an angle makes easy contact with the classical world, but 
once more considers the phase as a mere state parameter 
instead of a full quantum variable. 

The third major concept emphasizes the idea that 
phase is a physical property and, by any orthodox pic- 
ture of quantum mechanics, must be associated with a 
selfadjoint operator (or at least with a family of positive 
operator- valued measures). In this vein, phase opera- 
tors have been constructed via a polar decomposition for 
qubits and qutrits [Tof. 

The main goal of this paper is to look at the funda- 
mental problem of properly defining phase from quite a 
different perspective. On closer examination, one imme- 
diately discovers that the idea of complementarity is at 
the root of all the previous approaches: phase is comple- 
mentary to some amplitude, by which we loosely mean 
that the precise knowledge of one implies that all possi- 
ble outcomes of the other are equally probable Jul] . This 
idea of unbiasedness leads directly to introduce mutu- 
ally unbiased bases (MUBs) .12], which, for a variety of 
reasons, are becoming an important tool in quantum op- 
tics |13| . It is known that the maximum number of such 
bases cannot be greater than d + 1 and that this limit 
is certainly reached if d is a power of a prime [T^j. It 
is not known if there are nonprime-power values of d for 
which this bound is attained. We shall be not concerned 
with this problem in this paper, and assume that we are 
always working in a prime dimension. 

It is essential to recall that complementarity for the 
position-momentum pair is implemented by the Fourier 
transform, which exchanges both operators. Using the 
ideas introduced in Ref. [15j . we construct d+1 disjoint 
classes of maximally commuting unitary matrices (each 
set having d — 1 operators). We then note that one of 
these classes consists solely of diagonal operators (which 
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we can relate to inversions) that can be mapped, using 
the finite Fourier transform, to operators acting cyclically 
on basis states. 

This perspective leads to a natural notion of phases as 
complementary to inversions. The advantage of this ap- 
proach is that it does not rely on polar decompositions 
or semiclassical arguments and provides a clear under- 
standing of the behavior of such basic variables. 



II. MULTICOMPLEMENTARY OPERATORS 
FOR FINITE-DIMENSIONAL QUANTUM 
SYTEMS 

The objects we study in this paper are quantum sys- 
tems described in a <i-dimensional Hilbert space Tid- We 
recall 0] that two different orthonormal bases A and B 
are said to be mutually unbiased if a system prepared in 
any element of A (such as \a)) has a uniform probability 
distribution of being found in any element of B 
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(2.1) 



for all a E A and all b G B. As anticipated in the In- 
troduction, we shall be concerned solely in cases where 
d is a prime number, as then we know that there are 
d + 1 MUBs. For dimensions which are power of a prime 
the argument can be easily extended with some modifi- 
cations 0| 

If \n) [n = 0, . . . , d— 1) is the standard (computational) 
basis in Tid, we introduce the generalized Pauli matrices 
X and Z by the following action: 



where 



X\n) = \n+l), 
Z\n) = oj n \n), 

uj — exp(27ri/d). 



(2.2) 



(2.3) 



Note that throughout this paper addition and multi- 
plication must be understood mod d. These operators 
X and Z , which are generalizations of the Pauli ma- 
trices, were studied by Patera and Zassenhaus in 
a purely mathematical context, and have been used re- 
cently by many authors in a variety of applications |18| . 
Under multiplication, they generate a finite subgroup of 
SU(rf), known as the generalized Pauli group, and obey 
the finite-dimensional version of the Weyl form of the 
commutation relations: 



ZX = ujXZ. 



(2.4) 



It is easily shown that the eigenvectors of X and those 
of Z satisfy (JOJ. 

To simplify as much as possible the following compu- 
tation, we introduce the following labeling scheme: let 



Since we shall also need powers of these operators, we 
denote by <L% k the set 

£x k = {%k,%i, ■ ■ ■ , X^ 1 }. (2.6) 

The d—1 operators in the class £x k clearly commute one 
with another and therefore represent a maximal set of 
commuting operators. 

Following the ideas in Ref. consider now the fol- 
lowing set (each containing d + 1 operators) 



S — {Xo, Xi, . . . , Xd} 

= {z,x,xz.,...,xz d - 1 }. 



(2.7) 



By virtue of the relation 12. 4|) , any two operators in this 
set are complementary, in the sense that their eigenvec- 
tors satisfy the unbiasedness condition (|2.1|l . Further- 
more, a complete MUB is obtained by constructing every 
eigenvector of every element in S, so we refer to S as a 
maximal set of multicomplementary operators. 

In the case of the standard position-momentum com- 
plementary variables, their eigenvectors form bases re- 
lated by the Fourier transform. The finite-dimensional 
Fourier transform can be defined as 191 



F 



1 



d-l 



4= E ^ 

s/d ^ 



\n){n'\ 



n,n' — 

with the properties 

FF t = F^F = 1 



F 



1. 



(2.8) 



(2.9) 



Using this definition one can check that X and Z are 
indeed Fourier pairs 



X = F^ZF. 



(2.10) 



There exist also an operator V that transform X — > XZ k . 
It has a diagonal form: 



d-l 



,-(n-n)(d+l)/2 



n)(n\ 



so that 



XZ k = v j(k xv k 



(2.11) 



(2.12) 



where we have assumed an odd-prime dimension (the 
case d = 2 need minor modifications, as we shall see 
in next Section). 

In physical applications, only d—1 populations can 
vary independently in a d-level system. In consequence, 
it is usual to work with d—1 traceless operators hj that 
measure population inversions between the correspond- 
ing levels, i.e., 



hj — Sjj 



S 



(2.13) 



where 



Xq 



Xjt 



XZ 



fc-i 



(2.5) 



Mil- 



(2.14) 
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These hj (usually known as the Cartan-Weyl generators) 
constitute a maximal Abelian subalgebra. Note that the 
diagonal operators in the class £x = {Z k } are linear 
combinations of hj , so both can be used indistinctly. 

On physical grounds, we expect phases to be comple- 
mentary to inversions. But inversions are invariant under 
phase shifts: if 

U(<p) = exp (-iY^VjhA , (2.15) 

where ip denotes {tp\, . . . , ipd-i), then 

U\<p)h k U{cp) = h k . (2.16) 

Thus, we can construct a continuous family of operators, 
all complementary to inversions, by conjugating any op- 
erator in the class <Lx k by U((p). In particular, the diago- 
nal operator V, which maps X to XZ k as per Eq. (|2.12(l 
(and that is also of the form lj2.15|) for a definite choice 
of the parameters), allows us to construct any start- 
ing with Xi = X. For this reason, and without any loss 
of generality, we shall henceforth focus on the elements 
of the class C^i = {X, X 2 , . . . , X^ 1 } to discuss general 
properties of complementary phase operators. We thus 
define d — 1 families of operators representing the expo- 
nential of the phase by 

E k {p) = U^{p)X k U{p), fc = l,...,d-l, (2.17) 

obtained by successive powers. 

If \s) is an eigenstate of hj with eigenvalue hj s , then the 
expectation value of E k (cp) on an arbitrary state = 
J2 s c s\s) is simply 



(E k (p)) = E <+kCs U: +k (cp)U s (<p), 



where 



U s (p) = exp -i^2ipjhj s 



(2.18) 



(2.19) 



This allows us to introduce an operator kernel 

d-1 



r%>) = 



i 



(27T) 



d-1 



1 + J2( E * '(V))* xk 



fc=l 



(2.20) 



which is properly normalized and generates all the mo- 
ments through the relation 

(J \d-l 

i E (<P)) = ~ — 2 TrPM^'], (2.21) 

so the phase distribution is obtained as 



(2.22) 



We note in passing that any positive operator-valued 
measure of the general form 



Ate) 



1 



(2tt 



id-l 



d-1 



1 + E 7fe £ fc te) 



fc=l 



(2.23) 



which can be associated to ()2.22(l . satisfies the usual re- 
quirements of real valuedness, positivity and normaliza- 
tion and possesses the obvious property 

e^Atex, . . . , ? d i)eM>h = Ate x , . . . , iP j+ ( p' j , ^.i), 

(2.24) 

which meets the usual requirements of complementar- 
ity [Ij. 



III. APPLICATION: QUANTUM PHASE FOR 
FINITE-DIMENSIONAL SYSTEMS 



A. The case of qubits 

To fully appreciate the details of the method, we shall 
work out some relevant examples. First, we focus on the 
simplest case of a two-dimensional Hilbert space TI2, and 
a state space that coincides with the sphere S2. 

The basic operators are the standard Pauli matrices 



X = 2a x 



such that 



1 

1 



2ct 2 



1 

-1 



OVCr 2 = —O z O x 



(3.1) 



(3.2) 



The transformation a z — > cr x is accomplished by the finite 
Fourier transform 



F 



J_/l 1 

M 1 -1 



(3.3) 



However, it is impossible to find a unitary transformation 
V such that a x — > a x a z . For this reason, instead of a x a z 
the matrix a y = ia x a z is used, so that a y = V'a x V, 
where V is the unitary operator 



V = 



1 

-i 



(3.4) 



If qa describes a point in the Bloch sphere S% 
in an arbitrary direction parametrized by n = 
(cos if a sin $ a , sin <pa sin §a , cos *&a ) , let A be 

A = n A <r = R{$a, Pa) o z R-\d A ,VA) , (3.5) 

where 



R{"&, <p) = exp 



— (cos <p a x - sin ip a y ) 



(3.6) 
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The condition of complementarity between A and a 
generic operator B [expressed also as in Eq. (|3.5[) ] can 
be written as 1211 



n A ■ n B = ; 



(3.7) 



that is, the subspace spanned in S 2 by is orthogonal 
to that by ub- We have then a one-parameter set of 
complementary operators of the general form 



B = n B ■ cr , 



(3.8) 



where the unit vector satisfies <|3.7[1 . which is equiva- 
lent to 



cot$? 



tani?^ cos(ipB — 'Pa) ■ 



(3.9) 



In particular, the complementary set to the inversion a z 
consists in the one-parameter family 



E(ip) = cos ip g x — sin ip <j y — 



e iV 
e- zv 



(3.10) 



where ip represents a reference phase. This in fact agrees 
with the exponential of the phase operator obtained via 
a polar decomposition [Toj | . 

According to the approach developed in this paper, 
we have now one family of phase operators that can be 
constructed as 



E(tp) = esxp(iipcr x /2)(r x exp(—i(pa z /2). 



(3.11) 



This coincides also with (|3.10|) . and can be recast in the 
suggestive form 



E(<p) =e{<p) t Fa, 
where t denotes the transpose and 
1 



V2 



( e: " ] 






e 1 ' )' 







(3.12) 
(3.13) 



This result confirms in this simple case the complemen- 
tary character of E(tp) obtained via Fourier transform. 
For a pure state such as 

) ■ 

with < i3 < tt, 0<x< 2tt, the average value of E(<p) 
is 

(E(<p)) = sini9cos(x + tp). (3.15) 

The main features of this description are obviously inde- 
pendent of the reference phase tp. 

B. The case of qutrits 

For a three-dimensional Hilbert space H3 the basic op- 
erators are 



X 



















I 







h 





UJ 






1 











- 2 / 



(3.16) 



and uj = exp(27ri/3). We have four classes of disjoint 
traceless operators, each containing two commuting op- 
erators 



<£*„ = {Z,Z 2 }, £ Xl ={X,X 2 }, 



(3.17) 



<t Xl = {XZ, (IZ) 2 }, £ Xs = {XZ 2 , {XZ 2 ) 2 }. 
The discrete Fourier transform is 

111 



F = I 1 uj uj 2 
V3 \ 1 ^2 u 

and V is the diagonal unitary matrix 
V 



(3.18) 




(3.19) 



In this case we have two Cartan operators associated 
with the two independent inversions 



















• 


-1 


• 


hi = 


I 


1 


i 



















which can be easily expressed as linear combinations of Z 
and Z 2 . Note that X and X 2 correspond to two different 
physical situations: in the computational basis X acts as 
X\n) = \n+ 1), while X 2 \n) = \n + 2). 

Thus, we have two families of commuting phase oper- 
ators: 



E(pi,ip 2 ) — e i (^ lhl +^ 2h2 ^ X e^ i( - Vlhl+(p2h2 ^ 
E 2 (pi,<p 2 ) — e i ^ lhl+V2h2 ^ X 2 e - i (fit>-i+<P2h2) _ 



(3.21) 



In this particular case, they essentially coincide because 
E 2 (ipi,tp2) = E'{<px, (p 2 ). The phase operator E(<pi,<p 2 ) 
can be represented in a form similar to i|3.12[l . namely 



E(ipi,ip 2 ) = e((pi,ip 2 ) t F2L, 



(3.22) 



with 



e(<pi,<P2) = ■ 
Its explicit form is 



/ e H.V2-2<Pi) \ 






1 e i{<P\+<P2) J ) 


X = 


\xz 2 


\ e i(¥>i-2y 2 ) / 







(3.23) 



e *(<Pi+V2) 

E(<pi,ip 2 ) = ( e^ 2 ^^ 
e *(ipi-2¥>2) q 

(3.24) 

where again ipi and p 2 are reference phases. The average 
value of E(ip\, (p 2 ) in an arbitrary pure state such as 



cos(tf/2) 
sin(0/2) cos(£/2) e iXi | , 
sin(i?/2) sin(£/2)e 4X2 



(3.25) 
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is simply 

+ isinl9{sin(^/2)e i [ X2+ ^ 1+¥,2 )] 
+ cos{^/2)e-' l ^ 1+ ^ 1 ~^}. (3.26) 
We can observe that when only two levels are involved, 

(Efatfr)) = ^sinde-^+P^-^, (3.27) 

which measure the relative phase between the states 1 
and 2, and depends only on a single effective phase 
2(fi — (f2- Since E 2 (ip 1, 1P2) = E'((pi, (P2), which includes 
all the possible moments, the phase distribution (|2.22|) 
remains constant at the direction 2<p\ — ip 2 = constant. 

It is instructive to compare our construction of phase 
operators with the more common algorithm, based on 
the polar decomposition of Syz, S23, and S13, defined 
in Eq. I|2.14[) . The polar decomposition of Sji implicitly 
focuses of the i — > j transition, without analyzing the 
role of the third "spectator" state of the system. This 
observation can be used to explain the lack of uniqueness 
in the polar decomposition of singular operators like S12, 
for instance. This approach produces a phase operator 
E12 of the form 

/ 1 0\ 

Ei 2 = \ x y , \x\ 2 + \y\ 2 = 1, (3.28) 
V y* -x* J 

where S12 = -E12-R12 and R\ 2 = VS12S21 is the "mod- 
ulus" . The only constraints on _E 12 are imposed by the 
requirement of unitarity. 

A particular solution to Eq. I|3.28|l is obtained by 
"isolating" the third state from the first two by choos- 
ing y = and thus x = e lVl2 . Since (512, S21, hi) 



span a su(2) subalgebra, this choice amounts to limit- 
ing the action of the phase operator E12 to a specified 
su(2) subspace. A similar argument holds for E23, with 
(S23, S32, /12) spanning another su(2) subalgebra. 

This perspective in terms of polar decomposition and 
transitions leads to phase operators that fulfill the re- 
quirements of complementarity only between pairs of 
states involved in each transition. The corresponding 
positive operator-valued measure obtained in this con- 
struction is still of the general form found in Eq. I|2.23|l . 
However, this restricted point of vies is to be contrasted 
with the approach of this paper, where complementarity 
is imposed for the three-level system as a whole. 

IV. CONCLUDING REMARKS 

Mutually unbiased bases are a primitive of quantum 
theory, as they embody the importance of the super- 
position principle. In this paper we have used them to 
develop a comprehensive quantum theory of the phases 
as complementary to inversions in finite-dimensional sys- 
tems. 

The construction presented in this paper is devoid of 
any ambiguity associated with the non-uniqueness of po- 
lar decomposition of ladder operators. In prime dimen- 
sions, phase operators and inversions are elegantly re- 
lated by a finite Fourier transform, much like positions 
and momenta are related by an ordinary Fourier trans- 
form in infinite-dimensional systems, and provides an ap- 
pealing way of treating a concept as central as the phases 
of a system. 
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